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We consider the nonlinear propagation of electrostatic wave packets in an ultra-relativistic (UR) degen- 
erate dense electron-ion plasma, whose dynamics is governed by the nonlocal two-dimensional nonlinear 
Schrodinger-like equations. The coupled set of equations are then used to study the modulational instability 
(MI) of a uniform wave train to an infinitesimal perturbation of multi-dimensional form. The condition for 

1 /3 

the MI is obtained, and it is shown that the nondimensional parameter, /3 oc Xcn^ (where Ac is the reduced 
Compton wavelength and ng is the particle number density) , associated with the UR pressure of degenerate 
electrons, shifts the stable (unstable) regions at ~ 10"^° cm~^ to unstable (stable) ones at higher densities, 
i.e. no > 7 X 10^"^. It is also found that the higher the values of no, the lower is the growth rate of MI with 
cut-offs at lower wave numbers of modulation. Furthermore, the dynamical evolution of the wave packets 
is studied numerically. We show that either they disperse away or they blowup in a finite time, when the 
wave action is below or above the threshold. The results could be useful for understanding the properties 
of modulated wave packets and their multi-dimensional evolution in UR degenerate dense plasmas, such as 
those in the interior of white dwarfs and/or pre-Supernova stars. 
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I. INTRODUCTION 



As is known, compact astrophysical objects, e.g. massive white dwarfs and/or the core of pre-Supernova stars, 
are supported by the pressure of degenerate electrons, i.e., in their interiors the particle number density is extremely 
higlPj^. So, the Fermi energy can be much larger than the thermal energy (e.g. typical electron Fermi energy, 
Ep ~ 1 MeV corresponds to a Fermi temperature Tp ~ 10^'^ K). Thus, the thermal pressure of electrons may 
be negligible compared to the Fermi pressure. Such degenerate electrons may, however, be either nonrelativistic, 
somewhat relativistic or ultra-relativistic. In the latter case, the speed of electrons can approach the speed of light 
in vacuum (c), and the equation of state for electrons can then be written a^^Hll = (S/tt)^^'^ {hc/8)nt^^ whenever 
Ep ^ rrieC^ . Here, h is the reduced Planck's constant, is the electron number density and rUe is the electron mass. 
Such ultra-relativistic degenerate (UR D) el ectrons are, indeed, ubiquitous in many other astrophysical environments 
including neutron stars and magnetar^^'^'^. 

Now, in dealing with those compact objects like white dwarfs, their interiors may be considered as a plasma system 
(e.g., carbon-oxygen white dwarf in a thermonuclear Supernova explosion) consisting of positively charged ions (nuclei) 
providing almost all the mass (inertia) and none of the pressure, as well as electrons providing all the pressure (restoring 
force), but none of the mass (inertialess) . So, interiors of such compact objects provide us a cosmic laboratory for 
studying the properties of such plasmas as well as nonlinear collective oscillations under extreme conditions, i.e. at 
the relativistically degenerate dense states with higher densitie^ (10^ — lO^g/cm'^). Recent investigations along these 
lines indicate that such URD dense plasmas can support the propagation of solitary waves as well as double layers at 
different length scales of excitation^. 

On the other hand, the nonlinear propagation of wave packets in plasmas is generically subject to their amplitude 
modulation due to the carrier wave self-interaction, i.e., a slow variation of the wave envelope due to nonlinearities 
(see, e.g.l^ti^. Under certain conditions, the system's evolution shows a modulational instability (MI), leading to the 
formation of envelope solitons through the localization of wave energy. Such solitons, governed by a (l+l)-dimensional 
nonlinear Schrodinger equation (NLSE), are the result of a balance between the nonlinearity (self- focusing) and the 
wave group dispersion. When such a balance can be maintained dynamically, the solitons may exist even under 
strong perturbations. On the contrary, the multi-dimensional [(2-1-1) or (3-|-l)-dimensional] NLSE with cubic and/ 
or quadratic (nonlocal) nonlinearities may no longer be integrable. In this case, the self-focusing effect can dominate 
over the disper sion (beyond a threshold intensity) leading to the formation of wave collapse or blowup of the wave 
amplitudfPE2l. 

The important issue for the multi-dimensional NLSE may be that for a wide range of initial conditions, the system 
often exhibits collapse in which a singularity of the wave field is formed in a finite time, instead of a stable mode 
propagation. A collapsing wave packet thus self-focuses in shorter scales and with higher amplitudes until other 
physical effects intervene to arrest it. Such collapse ph enom enon plays an important role as an effective mechanism 
for the energy localization in various branches of physic J^^^. In this context, it could be a central problem of finding 
a proper initial condition which leads to wave collapse. This is, one of our goals of the present work in the description 
of Davey-Stewartson II (DS II)-like equation&i^j for the wave packets. Several authors have investigated such collapse 
theory, e.g. in the context of hydrodynamicJSl^ in nonlinear optic^^^, in Bose-Einstein condensateJ^^ as well as in 
plasma physicJ^. However, to our knowledge, there are few theoretical studies for the MI (see, e.^^^^^ and multi- 
dimensional evolution of electrostatic wave packets (see, e.g.l^ in plasmas, which are described by the DS Il-like 
equations. The latter generalize the (l+l)-dimensional NLSE with a nonlocal (quadratic) nonlinear term associated 
with the static field in the plasma. 

Furthermore, since quadratic nonlinearities are known to be collapse-free, multi-dimensional NLSE, where the cubic 
and quadratic (nonlocal) nonlinearities compete, can also support coherent structures, i.e. dromion-like solutions 
(unlike solitons the dromions can have inelastic collisions and can transfer mass or energy) which may either decay 
due to the dispersion to be enhanced by the static field or exhibit blowup due to nonlinearity, in a finite time^*^ 
However, their applications in physical systems, especially in plasmas are not yet fully understood or less developed 
till now. 

In this article, our purpose is to consider the propagation of two-dimensional (2D) electrostatic wave packets (EWPs) 
in a URD dense plasma. We provide a general criterion for the MI of a plane wave packet as well as the instability 
growth rate. We show that as one approaches the higher density regimes, the stability of the wave increases with 
lower growth rates at lower wave numbers of modulation. Furthermore, the 2D evolution of the nonlocal NLSEs 
exhibit dromion-like solutions which either decay by the wave dispersion or blowu p due to wave nonlinearity in a 
finite interval of time. The latter are in qualitative agreement with the result^^SllU already found in DS II equations 
for water waves with finite deptlP. 
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II. BASIC EQUATIONS AND DERIVATION OF THE EVOLUTION EQUATIONS 



Let us consider the propagation of EWPs in a 2D dense plasma composed of inertialess URD ultra-cold electrons 
and inertial ultra-cold ions. Any speed involved in the plasma flow is assumed to be much lower than the ion- 
acoustic speed. At equilibrium, both species have equal number density, say Hq. Assuming further that the plasma is 
coUisionless and unmagnetized, the basic normalized equations then reacP 

^+V.(n,v)=0, (1) 

|^ + (v.V)v = -V</), (2) 

= V<t>~^Jnt/\ (3) 
V^^ = He - n„ (4) 

where V = {d/dx,d/dy), Ui is the ion number density normalized by riQi v = {vxTVy) is the ion velocity normalized 
by Ci — ^ ZiVfif^c? jnii with Zi denoting the ion charge state, rrii is the ion mass and c is the speed of light in vacuum. 
Also, is the electrostatic wave potential normalized by ZiirieC^ /e with e denoting the elementary charge. The space 
and time variables are respectively normalized by the screening length = yJm^Jz^^jAimt^ and the ion plasma 
period, u)^,^ — y^mi/ATmoZfe^ . Moreover, in Eq. /? = Ac\/?io/727r where Ac = h/mc is the reduced Compton 
wavelength, and we have used the same equation of state for Pg as described in the previous section relevant for URD 
electrons in dense plasmas. Furthermore, the nondimensional parameter, /3 ^ 1 for uq ^ ric ~ 3.94 x lO'^'^cm"'^ and 
/? = 1 for uq ~ fic, i.e. higher values of j3 represent the higher-density regimes. 

Next, we consider the time evolution of a wave packet of electrostatic perturbations which occur along the a;-axis. 
The initial wave packets of perturbation are basically modulated by the nonlinear carrier wave self-interactions. Then 
if one observes the wave packets from a coordinate frame moving with a group speed Wg, to be obtained from the 
linear dispersion relation (as given below) of Eqs. Q-Q, then the time variation of the wave packets looks slow, and 
so the space and the time variables can be stretched as 

£, ^ e{x - Vgt),ri ^ ey,T ^ e^t, (5) 

where e is a small parameter representing the strength of the wave amplitude. We are interested in the modulation of 
a plane wave as the carrier wave with wave number and frequency k and uj respectively, then following the standard 
reductive perturbation technique (RPT) (see, e.g.^i^, the dynamical variables can be expanded as 



oo oo 



ne,^ = l + ^e" n["^^),(^,?7,T)exp[i(fca;-wt)Z], (6) 

n— 1 /— — oo 

oo oo 

Vx.y = XI ^" ^) [*(^^ ~ ^*)'] ' 

n— 1 oo 

oo oo 



n=l /=— oo 



where "^("y)/ ^^'^ should satisfy 5*1" = sj"^* because of the reality condition for the physical variables. 

Here, the asterisk denotes the complex conjugate of the corresponding quantity. Note that the group speed, Vg is now 
normalized by Ci (smaller than the ion-sound speed), the wave frequency w and the wave number k are normalized 
by ujpi and A^ respectively. 

We now substitute the expressions (|6|-(|8]) into the basic Eqs. Q-Q and equate the terms in different powers of e. 
We shall, however, omit the detail calculations, since the procedure is quite standard and follows the usual RPT. In 
the lowest order of e, we obtain for n = I, I = 1 the linear dispersion relation [since we are considering the modulation 

of a plane wave, n^(l.\\p v^,^^ ^ and (Z)!^"* are all set to zero except for I — ±1] 



2 /^^ fn\ 

We find that the dispersion equation ^ has the similar form with that derived by Kako et aP^l for classical plasmas 
comprising cold ions and isothermal electrons whenever one replaces /3 by unity. However, one should note that the 
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case /? = 1, which represents the higher density regimes where Uq = ric, is not apphcable for classical isothermal 
plasmas. Equation ^ shows that the wave always propagates with a frequency below the ion-plasma frequency 
regardless of the values of /? and k. In the short- wavelength limit, the frequency of the wave approaches unity (the 
upper limit of the normalized lo), whereas in the long-wavelength limit, it approaches a zero value. Furthermore, 
the limit /Jfc^ ^ 1 is not admissible because otherwise, oj will be small to provide weak dispersion, and the soliton 
formation (Korteweg-de Vries soliton) of the carrier wave will be a lower-order process than the MI of the envelope 
to be studied here. Moreover, since k is normalized by the inverse of the screening length A^, the values of fc > 27r are 
also inadmissible, otherwise the wavelength would become smaller than the screening length. As a result, the plasma 
collective behaviors might disappear. 

Now, proceeding in the same way as of Refs.'^^'^, i.e., considering the second harmonic modes obtained in terms of 

(2) fl) (1) ' (1)' 2 

(p\ , d(j)\ /d^, d(p\ /dr] for n = 2, Z = 1 and (p\ for n = 2, Z = 2 as well as the zeroth-harmonic modes appearing 

due to the nonlinear self-interaction of the modulated carrier waves, and finally considering the equations for n = 3, 
I = 1, we obtain the following 2D nonlocal NLSEs for the propagation of modulated wave packets in URD dense 
plasmas 



g^2 Q^2 g^2 



-P2 + Qi|0|'<^ + Q2V'0 = O, 



(10) 

(11) 



where (j) = (t)i , = / d^Vy^drj and the coefficients written only in terms of k and the parameter /3 are given by 



Pi - 

Qi = 



Q2 

S = 



3k/3 



3/2 



2(l + /3fc2)^/2' 



ky/lTW 

k{l + /3k^) 
k{5~ 3/3fc2) 



3/2 



fp/2 



3 + (3k^ 



P2 



' 2fc(l + /3fc2)3/^' 
7^^ 



(12) 



6V/?(l+/3fc2)3/2' 

1 + I3k^ 



3(l-l-/3fc2)2 

\ , /3-M5-2//3^) 
3(l + /3fc2)3 

R ^ 1 



2 + 3k 



(3^k^^ 
fc(-2 + 3^ + 3/32fc2) 

9/35/2 (1 ^ /Jfc2)3/2 



(l + /3fc2)'^ 



3il + l3k^f 



(13) 
(14) 
(15) 



We note that since, Pi < 0, P2 > and R> 0, the NLSEs ( |10[ ) and (11) are in the form of second type (hyperbolic- 
elliptic) Davey-Stewartson, i.e. DS Il-like equations. Furthermore, 5* is always negative, Q2 ^ for /3fc2 ^5/3, 
Q2 = for f3k'^ = 5/3 and Qi can be either positive or negative depending on the values of both k and /3 . In the 



case Q2 = for which Eqs. (10) and (11) are reduced to 2D NLSE, the modulated wave can be shown to be stable. 



This will be clear from the next section dealing with the MI. We will see that for values of /3 > 1 (which corresponds 
to higher density regimes) and small wave numbers, i.e. fc < 1 fulfilling /3fc2 < 3, the wave tends to become stable. 

Thus, we have obtained a new set of nonlocal NLSEs, which describe the slow modulation of EWPs in 2D URD 
dense plasmas. The coefficients Pi, P2 appear due to the wave group dispersion and the 2D evolution of the EWPs. 
One of the nonlinear coefficients, Qi (cubic) is due to the carrier wave self-interaction originating from the zeroth 
harmonic modes (or slow modes) and the other nonlocal coefficient (quadratic) Q2 comes from the coupling between 
the dynamical field associated with the first harmonic (with a 'cascaded' effect from the second harmonic) and a static 
field generated due to the mean motion (zeroth harmonic) in plasmas. 



III. MODULATIONAL INSTABILITY AND GROWTH RATE 



In this section, we consider the modulation of a plane wave solution of Eqs. (10) and (11) for 
amplitude 4>o- The boundary conditions, namely (/>, -0 — >■ as ^, 77 — >■ 00 used before must now be relaxed 



with a constant 
since the wave 



packet is still not modulated and the solution is not unique. Note that the choice of '0 = ■00 is immaterial as the stability 
criterion does not depend on it. Thus, we can represent the plane wave solution as = 0o exp [i{k^^ + fc^ry — Ar)], 
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and ip = i^o with A 



Pikl + P2^ 



Q2'(po, where k^, fc^ 



tpo are all real constants. Next, to study 
the stability of this solution we modulate the amplitude as = (0o + </'m) exp [i{k^£, + kr/r] — At)], ip = ipo + i^m, 
where ■(/'m oc Re [expi(iiri^ + ?7 — ^t)] and K, Q. are respectively the wave number and the wave frequency of 
modulation. Looking for the nonzero solution of the small amplitude perturbations, we obtain from Eqs. ( 10 ) and 



(111 the following dispersion relation for the modulated wave packet 



where 



K, 



= {PiKl + P^Kl) (1 - Kl/K'') , 



{Pi+P2e){R+e) 



(16) 



(17) 



in which C = tan0 — K2/K1. Equation (16 1 shows that the MI sets in for a wave number satisfying K < Kc or, for 

(18) 



all wavelengths above the threshold, Ac — 2t:/Kc, provided the right-hand side of Eq. (17) is positive, i.e. 

{R + C^)Qi+SQ2 



{Pi+P2e){R+e) 



> 0. 



Since i? > 0, the condition ( 18 ) reduces to 



A=[{R + e) Qi + SQ2] {Pi + P2C) > 0. 



(19) 



The wave packet is otherwise (i.e., for K > Kc) said to be stable under the modulation. The instability growth rate 
(letting f2 = iT) can be obtained from Eq. ( 16 1 as 



(Pi + P2C' 

i + C" 



1. 



(20) 



Clearly, the maximum growth rate, achieved ai K = is Fn^ax — \{R + C^) Qi + SQ2\ I'/'oP/ {R + C^)- 

Next, we numerically investigate the condition of MI given by Eqs. (18) or (19) as well as the instability growth 
rate, T given above. The condition of MI not only depends on the carrier wave number k and the density dependent 
parameter (3 arising due to the ultra-relativistic pressure of degenerate electrons, but also on the obliqueness parameter 
of the modulational wave number K with the ^-axis due to 2D perturbation. The stable (A < 0) and unstable (A > 0) 
regions are thus shown in the kO plane in Fig. 1 for different values of j3 that correspond to different density regimes. 
We find that the stable and unstable regions are completely divided into two parts at comparatively lower as well as 
higher density plasmas, i.e. at no lO"^" cm^^ [see Fig. 1(a)] and at no = 5 x 10'^'^,7 x lO'^'^ and 10^'' cm~'^ [see 
Figs. l(d)-l(f)]. Separation of such regions are also observed at the intermediate densities [Figs. 1(b) and (c)]. As 
seen from these figures, there arc basically four regions (two for each) for the stable and unstable waves, and for a 
wide range of values of k and tan 6'. From Figs. 1(a) and 1(e) or 1(f), one observes that some part of the regions in 
which the wave was stable (unstable) at lower density now shifts to unstable (stable) region at higher values of the 
same. Obviously, the parameter responsible for this shift is j3 due to the consideration of ultra-relativistic degenerate 
electrons. 

Again, P2, R are always positive and Q2 > ^ for /3fc^ < 5/3. Also, 5* < and P2 >| Pi | for pk'^ > 1/2. Then, in the 
regime satisfying 1/2 < < 5/3, the instability condition (18) or (19 1 depends mainly on the sign and magnitude 
of the coefficient Qi. Physically, this implies that the EWPs become stable or unstable in the said regime due to the 
nonlinear self-interactions originating from the second harmonic modes as well as from the zeroth harmonic modes 
(or slow modes). In particular, under the horizontal modulation {Ki = 0) the instability condition ( jlS] ) reduces to 
A = Q1/P2 which can be either positive or negative depending on the sign of Qi (since P2 > 0). We also find that 
for a particle density to vary in 10'^^ — 10^'^ cm~'^ and for carrier wave numbers fc < 1, Qi is negative, and then 
the modulated wave is said to be stable. On the other hand, under the longitudinal modulation, i.e. K2 — 0, we 
have A = {QiR + SQ2) /Pi, which may, however, be positive almost everywhere in the above regime, and hence the 
instability. These results are in qualitative agreement with the classical ones considered before by Nishinari et ap3 in 
the description of DS II equations for EWPs. However, the general situation is quite different as clear from Fig. 1. 

Now, the MI growth rate can be calculated from Eq. (20 1 for a fixed value of the carrier wave number. Figure 2 



shows that higher the density regimes, the lower is the growth rate of instability with cut-offs at lower wave numbers 
of modulation. It seems that the parameter /3 plays almost the similar role of dispersion as the quantum diffraction 
associated with the Bohm potential plays for modulated wave packets in quantum plasmaJ^^^^^. This is expected as 
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(a) (6) (c| 




FIG. 1. (Color online) Contour plots of A =const. against k and tan 6'. The stable (shaded and/or wherever mentioned) and 
unstable (white and/or wherever mentioned) regions are shown in the feS-plane. The panels (a) to (f) represent respectively 
the regions corresponding to the densities, no (cm"^)= 10^°, 10^^, 10^^, 5 x 10^^, 7 x 10^^ and 10^*. The critical value at 
which the stable (unstable) regions are shifted to unstable (stable) ones is no — 7 x 10^"^ cm~^, and above which the stability 
region increases slightly with no [see Fig. 1(f)]. The panels (a), (e) and (f) show that the regions of stability and instability 
are divided into two parts at both lower and higher densities. 



we have seen from Fig. 1 that the waves tend to become stable at higher number densities. Thus, from this section 
we conclude that the number densities of ultra-relativistic ultra-cold electrons significantly modify the stability and 
instability regions in the fc^?-plane as well as reduce the instability growth rate for the wave packets. In the next section, 
we will numerically investigate the dynamical evolution of the EWPs that undergo MI in the regimes discussed above. 
We will see that the situation is quite distinctive from the one-dimensional case in which an exact balance between 
the self-focusing and the dispersion can happen to form envelope solitons on stable wave propagation. 



IV. 2D EVOLUTION OF THE NONLOCAL EQUATIONS 



In the evolution of EWPs described by the Eqs. (10 1 and (11) we first obtain some analytic conditions for the 



wave collapse to occur within a finite time. Note that since the evolution equations are of the DS Il-type, we can not 
derive the criteria that are sufficient to ensure collapse by the Virial theorem or els^2]_ However, we will present some 
conditions according to Berkshire and Gibbon'^. To that end, we first see that the integrals of motion are the wave 
action N and the Hamiltonian H where 



N = 



(21) 



n = 



Pi 



(^2 

2S 



R 



P2 



drj 



2S 



d^dr] 



d£,d-q, 



(22) 
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=6xl0^^cm"^ 



0.2 0.4 0.6 0.8 
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FIG. 2. (Color online) The modulational instability growth rate, F is shown with respect to the wave number of modulation, 
K for a fixed k = 0.6, ( — 0.2 and 00 = 0.1. This shows that the higher the particle density, the lower is the growth rate, F 
with cut-offs at lower K. 



where d^u/d^^ = ip. If / is the moment of inertia of a locahzed wave form, then by the Virial theorem we havtP 



Pi P2 



which upon integration gives 



(U 
d^ 



= 8nT + Ci, 



or, 



N N 

The constants Ciand C2 are then given by 

_/dI\ _ dl{0) 



, C2 = (!) \r=0 ^ 1(0) ■ 



(23) 



(24) 



(25) 



(26) 



> 0, Ci (0) < 
= 0, i.e. dl/dr = 



Thus, according to Ref.^^, the conditions for the collapse to occur in a finite time are (d^I/dr' 
0, C2 (0) > and that the angular momentum integral vanishes, J {t — 0) = 0. Furthermore, if Ci 
8Ht, then if > an initial waveform stretched along ^-axis will evolve into a structure stretched along the 77-axis 
for sufficiently large times. On the other hand, for H < 0, the wave will be stretched along the ^-axis. The detailed 
analysis is, however, beyond the scope of the present study, instead we will focus on the numerical evolution of the 
nonlocal NLSEs (fTol) and (II). We find that whatever may be the initial waveform (e.g., 2D Gaussian or close to the 



exact solution), the wave amplitude either decays due to dispersion or blows-up in a finite time due to nonlinearity, 
when the wave action to the initial condition is below or above the threshold. Furthermore, for a particular value of 
the power of the initial condition close to the exact solution, we find that the coherent structures, such as lumps, can 
also propagate almost without any change. 

The behaviors of the solutions of the DS II equations are quite well-known and have been investigated by many 
authors in the context of water wave propagation with finite depth (see, e.g. Rcfs.^'^^^). However, the solution is 
still not well understood or almost unknown in the context of plasma physics. To our knowledge, few authors have 
investigated the MI (see, e.g. Refs.'^^'^^ of a plane wave packet as well as some analytic solution (see, e.g.'^ of the 
DS Il-like equations in plasmas. 

In our numerical scheme we consider a space domain as [—15, 15] x [—15 x 15] with 150 grid points in every direction 
and time step St = 10~^. The parameter values which satisfy the MI condition (18) are taken as Pi = —0.317, 



7 




8 



(a) (b) 




FIG. 5. (Color online) Double focusing solutions that blowup in a finite time. The initial wave form is Gaussian with amplitude 
A = 8. (a) Initial state: r = 0, (b) Intermediate state: r = 0.35, (c) Before blowup: r = 0.356 and (d) Shortly after blowup 
(1(^1 ~ 10^): r = 0.358. 



P2 = 0.9327, Qi = -0.1686, Q2 = 0.2876, R = 0.3028 and S = -4.2086 corresponding to fc = 0.4 and no = 2 x 10^3 
cm^'^, i.e. f3 = 0.7985. For the sake of simplicity, we consider a symmetric Gaussian profile as the initial condition, 
i.e. (p — Aexp (^—^'^ /\Pi\ — ri'^ / P2) ■ We have tested our numerical results with other forms of initial conditions, 
namely (j) ^ A/ (^1 + S^'^ / \ Pi \ +Tf /P2) exp [i^"^ / Px +iTf /P2), which is a localized lump with algebraic decay and also 
(j) = [A/ (1 + ^^/iPil + ?7^/^2)] exp (2i77/v^), which approximates an exact solution of the DS 11-like equationJ^HHI]^ 
but observe the similar qualitative features as presented here. 



Figure 3 shows that for A — \^ the initial wave forms decay with time, and go to zero after a finite time. The 
corresponding contour plots are shown in Fig. 4, but for different times. In this case, the presence of the nonlocal 
static field enhances the dispersion of the profile and since the quadratic nonlinearities (nonlocal term) are known to 
be collapse free, the initial wave action N is not up to the mark for which the cubic nonlinearity dominates. However, 
since the effect of the static field is to modify the nonlinearity originating from the zeroth harmonic modes (or slow 
modes) , we expect that the coupling of the field ■0 and the wave field amplitude (j) can drastically affect the evolution 
dynamics for blowup. To elucidate it we have considered A = which is above the critical power N of blowup for 
a Gaussian pulse whose evolution is described by the 2D NLSE. Thus, Figs. 5 and 6 (contour plot) show that the 
initial waveform is not preserved, and wave singularity is formed at two points within a finite interval of time leading 
to double focusing effect. The wave packet thus blows-up in shorter scales with higher amplitudes. In this case, the 
self-focusing effects dominate over the dispersion as required. Evolution of such wave collapse could be an effective 
mechanism for energy localization in URD dense plasmas. 



A slightly different perspective is observed when considering the initial profile as 

(j) = [A/ (1 + I -Pi I +V^/P2)] exp (2i77/v^), i.e. close to the exact solution of DS 11 equationpl2ll. Figure 7 
shows that for A = 1 the coherent structure propagates along the ^-axis almost without any change (after a certain 
interval) of its amplitude. We can, however, precisely deal with such cases by choosing the initial condition as close 
to the exact solution of the NLSEs (10) and (11) as well as the proper wave action. 
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FIG. 6. (Color online) Contour plot {\<l)\ =const.) of a double focusing solution that blows up (as in Fig. 5) is shown at an 
intermediate time r = 0.35. 

V. CONCLUSION 

We have investigated the multi-dimensional modulation of an electrostatic wave packet propagating in an ultra- 
relativistic ultra-cold degenerate dense plasma. The dynamics of such wave packets is described by a coupled set of 
nonlocal nonlinear Schrodinger-like equations that involve a nonlocal (quadratic) nonlinear term. The latter appears 
due to the static wave field originating from the mean motion (zeroth harmonic) in the plasma. The equations are 
then used to obtain the instability condition for the modulation of a plane wave packet. It is shown that the density 
dependent parameter /?, which arises due to the ultra-relativistic pressure of degenerate electrons, shifts the stable 
(unstable) region at lower density {uq ^ lO"^" cm^'^) to unstable (stable) regions at comparatively higher densities 
(no ~ 10'^^ cm^^), and also that the stable and unstable regions are completely separated in such regimes. The 
latter can be achievable, e.g. in the interior of massive white dwarfs and neutron stars^. Furthermore, the instability 
growth rate is obtained and found to be lowered at higher number densities with cut-offs at lower wave numbers of 
modulation. This implies that the role of /3 is analogous with that due to quantum dispersion associated with the 
Bohm de Broglie potential for the MI of wave envelopes in quantum plasmaa^^iSl, 

We have also shown that the presence of the static field can drastically change the dynamical evolution of the 
wave packets quite distinctive from the one-dimensional wave packets or 2D NLSE. We found that when the initial 
condition is either a Gaussian pulse or a localized lump with algebraic decay, the wave amplitudes either disperse 
away to zero or blowup to infinity at singular points when the wave action (power) is below or above a threshold 
value. Such a blowup mechanism could be important for the energy localization in dense plasmas where the electrons 
are ultra-relativistic ally d egenerate. On the contrary, when an initial waveform is close to the exact solution of the 
DS Il-like equation I, the coherent structure propagates without almost any change. However, confirmation of 
this behavior needs further detail numerical investigation. 
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